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I.  INTRODUCTION 


Roll-lift  dispersion  of  spinning  missiles  is  caused  by  nonaveraging  of 

the  lift  vector  in  the  presence  of  moment  disturbances  that  perturb  the 

angle  of  attack.^  Lift  nonaveraging  can  occur  either  from  angle-of-attack 

variations  at  constant  lift  precession  rate  or  from  roll  rate  variations  for 

a  trimmed  missile  with  constant  lift.4  Control  of  lift  nonaveraging  should 

therefore  be  possible  by  controlling  either  the  lift  variations,  the  lift 

vector  precession  rate,  or  a  combination  of  both.  The  latter  method  was 

investigated  to  control  the  dispersion  of  an  untrimmed  missile  in  which  both 

angle-of-attack  and  lift  precession  rate  feedback  were  used  in  the  control 
3 

law.  In  another  investigation  to  control  dispersion  of  a  rolling  trimmed 
missile,  angle  of  attack  and  angle-of-attack  rate  feedbacks  were  applied  to 
pitch  and  yaw  control  moments  at  constant  roll  rate.^ 

In  this  report  we  consider  only  roll  rate  modulation  about  a  steady  roll 
rate  to  control  lift  nonaveraging.  Roll  modulation  complicates  the  analysis 
because  the  missile  equations  for  complex  angle-of-attack  motion  contain  roll 
rate  terms  in  the  coefficients.  The  equations  are  linear  for  constant  roll 
rate  but  become  nonlinear  when  the  roll  rate  is  time  dependent.  The  magnitude 
of  the  roll  coupling  depends  on  the  coordinate  system  used  to  describe  the 
complex  angle-of-attack  motion.  Both  body-fixed  and  aeroballistic  coordinates 
are  used  to  investigate  different  feedback  controls.  The  advantages  of  the 
different  coordinate  systems  are  described. 
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II.  ANALYSIS 


A.  FORMULATION 

The  equations  of  motion  of  a  spinning  missile  can  be  expressed  in 
different  coordinate  systems.  Shown  herein  are  the  differential  equations  for 
the  complex  angle-of-attack  motion  in  body  fixed  coordinates,  Eq.  (1),  and 
aeroballistic  coordinates,  Eq.  (2), 

6  +  [u  +  ip  (2  -  u) ]  6  +  [w2  -  p2(  1  -  u) 

+  ifp  +  p(v  -  v  ))U  «  im.  (1) 

5  +  [u  -  iyp]C  +  [ui^  -  ipv  ]  =  im  e1  '  pdt  (2) 

m  t 

for  small  angle  motion.  The  complex  angles  of  attack  in  these  two  coordinate 
systems  are  related  by 


5  -  «  e1®  (3) 

where  ♦  is  the  roll  angle  in  inertial  space,  defined  by 

♦  -  J  pdt  (4) 

and  mt  represents  pitch  and  yaw  moment  disturbances  about  the  y  and  z  axes 
divided  by  the  pitch  or  yaw  moment  of  inertia, 


M  +  iMz 

on  *  m  +  im  -  — ^ — -  (5) 

t  y  z  I  v  ' 

Disturbance  moments  can  occur  from  various  sources.  During  the  reentry 
phase  of  a  ballistic  missile,  asymmetries  occur  from  ablation  shape  change, 
which  can  cause  lift  variations  leading  to  dispersion  errors.  The  lift  forces 


do  not  average  out  for  each  roll  revolution  giving  rise  to  a  net  force  per¬ 
pendicular  to  the  flight  path.  The  sideways  motion  produced  is  called  the 
transverse  velocity  and  can  be  calculated  by  integrating  the  missile 
acceleration  in  the  cross  plane.  By  assuming  the  missile  to  be  a  point  mass 
in  this  plane  under  the  action  of  a  rotating  lift  force,  the  transverse 
velocity,  V  =  v  +  iw,  is  given  by 

V  =  V(0)  -  —  /  6ei-^pdtdt  (6) 

m 

o 


Lfl  C 

v  =  v(o)  -  ~r  !  yt  (7) 

m  ' 
o 

where  Eq.  (6)  is  in  body  fixed  coordinates  and  Eq.  (7)  is  in  aeroballistic 
coordinates.  The  lift  force  derivative  La  is  in  general  slowly  varying  and 
assumed  constant  over  the  time  duration  of  the  disturbances. 


The  nonaveraging  lift  effects  can  be  compensated  for  by  implementing  roll 
control.  The  roll  acceleration  is  defined  by 


l  +  i 
x  c 


(8) 


where  l  is  the  roll  disturbance  moment  and  l  is  the  control  moment. 

X  c 

The  transverse  velocity  is  a  measurement  of  missile  dispersion  and  is 
used  as  the  cost  function.  There  are  two  additional  quantities  which  are 
desirable  to  minimize  the  magnitude  of  the  applied  control  moment  and  the 
deviation  of  the  roll  rate  from  a  specific  value.  The  complete  cost  function 
is 


I  +  |v 


steady 

state 


I  +  w. 


t  - 

/  *cdt 


+  W 


Ap  - 

o 


P0)2dt 


(9) 


where  Wc  and  Wp  are  suitable  weighting  functions 


The  second-order  differential  equation  for  the  complex  angle  of  attack., 
Eqs.  (1)  or  (2)  involves  coefficients  that  are  functions  of  the  roll  rate  p. 

A  nonconstant  roll  rate  gives  rise  to  nonlinear  terms. 

B.  STATE  VARIABLE  FEEDBACK 

In  this  approach,  the  control  moment  SL  is  assumed  to  be  a  linear 
function  of  the  state  variables  of  the  form 

i 

~  =  AAa  +  BAB  +  Ca  +  DB  +  E/a  +  F/S  (10) 

x 

The  equations  of  motion  are  linearized  in  the  body  fixed  coordinates  by  taking 
small  perturbations  <5+,  p+  about  quasisteady  values  6q,  pQ 

A  =  6  +  6.  (11) 

o  + 

P  =  PQ  +  P+  (12) 

and  neglecting  higher  order  terras. 

The  resulting  linear  control  equations  in  terms  of  the  perturbations  are 
5+  +  [v  +  ipQ(2  -  p)]6+  +  [ u>2  -  P02(l  ~  p) 

+  ipo(v  "  VK  +  [2pop+(1  "  p) 

+  1(pj.  +  Pj.(v  '  v  ))]5  =  im^  (13) 

+  +  m  o  t 

P  =  Aa^  +  BB  +  Ci  +  DB  +  Efa  +  F fa  (14) 


Roll  moment  disturbances  are  assumed  zero.  The  transverse  velocity  is 
expressed  as 


(15) 


where  we  have  defined  X 

T 

X  *  /  £-  dt  (16) 

o  po 

and  changed  to  a  nondimens ional  roll  angle  variable  t  defined  by 

dr  *  pQdt  (17) 

We  assume  a  steady-state  condition  defined  by  the  initial  conditions  p  =  pQ 
and  6  “  6^  at  t  *  0.  Assuming  y  *  v  =  0  and  neglecting  smaller  terms,  the 
equation  for  the  angle  of  attack  is,  with  respect  to  t,  where  d/dt  *  pQ  d/dt, 


2  d2(6+)  d(6+) 

p  - ~ —  +  2  ip  — j -  +  oj  <S  =  im. 

o  dT2  o  dr  +  t 


(18) 


The  Laplace  transformation  of  Eqs.  (14)  and  (18)  results  in 


6+(s)  «  g(s)  +  io(s) 


(imT)  [(s2  +  u2/po2)  -  2  is] 
P02[  (®2  +  “2/P02)2  +  4  s2] 


(19) 


2  2 

psA+p  sC+E 
p(s)  -  — - - -  a(s) 

po  8 


2  2 

p  sB  +  p  s  D  +  F 

+  — - y~2 -  0(s) 

Po  8 


(20) 


Also  in  the  Laplace  transform  domain  is 


The  transverse  velocity  can  be  evaluated  in  terras  of  d(s)  and  p(s)  by  the 
following  procedure:  The  upper  limit  of  the  integral  in  Eq.  (15)  is  taken  to 
be  sufficiently  large  to  include  the  time  duration  of  the  perturbations. 
Without  loss  of  generality  this  limit  can  be  taken  as  »,  and  Eq.  (15)  can  be 
written  in  the  form 


Lft  00 

V  -  V(0) - —  /  F(t)  e1Tdt 

p  m  J 
o  o 


(22) 


where 


F(  t)  »  (So  +  d+)  e 


iX(  t) 


(23; 


From  the  definition  of  the  Laplace  transform 


F(s)  -5?  [F(t)]  h  J  F(t)  e~STdt 
o 


(24) 


the  integral  in  Eq.  (22)  is  observed  to  be  the  Laplace  transform  of  F(t)  with 
s  =  -i,  or 


v  *  v(0)  -  (F(s) ] 

p  m  S”— 1 

ro 


(25) 


where  F(s)  is  obtained  in  terras  of  d(s)  and  X(s).  By  expanding  the 
exponential  in  Eq.  (23),  F(x)  can  be  approximated  by 


F(t)  »  6  +  5 . ( x)  +  i6  (t) 

o  +  o 

<5  X2(  t) 

+  i5+(T)\(x)  -  -2— -  +  ...  (26) 

which  has  the  Laplace  transform  with  respect  to  t 

5 

F(s)  -  —  +  6(s)  +  id  X(s)  +  i 3! {d(x)X(x)} 
s  o 

-  -la?  {XZ(T)>  +  ... 


(27 


Eq.  (22)  can  now  be  written  as 


V 


V(0)  - 


p  m 
o 


fi(s)  +  ifi  X(s)  +  i.2  (5(t)X(t)} 
o 


-  -f  i^Cx^C t)>  + 


(28) 


The  velocity  caused  by  the  initial  angle  of  attack,  averages  out  for  each 
cycle  and  is  expressed  by 


-L.  6 

[— ] 

L  sJs=*-i 


p  m 

*o 


The  velocity  which  causes  dispersion  is  then  the  increment 


-L, 


AV  =  — -  [  6(  s)  +  ifi  X(s)  +  i  &  { 5(  t) X(  t)  } 

n  m  ^  n 


p  m 
*o 


-  ^{X2(t)}  + 


(29) 


(30) 


1 .  SOLUTION  FOR  CONTROL 


In  order  to  compensate  for  the  lift  force  asymmetry  caused  by  pitch  and 
yaw  moment  disturbances,  the  roll  rate  is  adjusted  throughout  the  cycle.  For 
example,  the  roll  rate  is  decreased  when  the  lift  force  is  in  the  direction 
opposing  that  of  the  dispersion.  The  necessary  roll  motion  is  found  from  Eq. 
(30)  by  setting  AV  equal  to  zero 


where 


AV  - 


[F(s)]s 


(31) 


F(s)  -  5(s)  +  16  X(s)  +  i  SB  { 5(  t) X(  t) } 
o 


--fa?  (x2(t)}  +  ... 
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The  solution  is  the  feedbacks  that  cause  the  real  and  imaginary  components  of 
F(s)  to  vanish  when  s  ■  -i 


F(-i)  =*  0 


Examined  first  is  the  case  where  there  is  an  initial  nonzero  angle  of 

attack,  6  .  This  is  caused  by  an  existing  asymmetry,  either  built-in  or 
o  ‘ 

previously  caused  during  flight.  The  lift  vector  produced  from  the  initial 
angle  of  attack  is  of  constant  magnitude  and  precesses  in  space  at  the  roll 
rate  frequency.  With  a  constant  body-fixed  angle  of  attack,  the  lift  force 
averages  and  the  missile  executes  a  spiral  motion  about  the  mean  flight  path. 

A  first-order  approximation  to  F(s)  is 


Fj(s)  =  <$(s)  +  i6QX(s) 


The  complete  expression  to  determine  the  gains  for  an  arbitrary  moment 
disturbance,  mt,  is  obtained  by  substituting  Eqs.  (19)  and  (21)  into 
Eq.  (33)  and  setting  F^-i)  ■  0  which  yields 


TT  -  [M]  { tG]  +  i [H] } [K] 


where 


[Ml  -  lmy  mz] 


(Cl  -  1°  ♦  1) - 

po°((o  -  1)^  -  4} 


po<o-l> 


0  -2po4  0  2 


(o  +  l) 


P03((c  -  l)2  -  4}  I  Q 


-Po^n  -2P^ 


0  -2p  -po  ( o—l ) 


0  2  0 


P0  (q-1)  0  -(o-l) 


IK]”  -  (A  B  C  D  E  FJ 


a 


The  magnitude  of  the  control  moment  must  be  minimized.  Examination  of 
the  solution  reveals  that  the  gains  are  inversely  proportional  to  the  initial 
angle  of  attack.  The  lift  force  resulting  from  the  angle  of  attack  provides 
the  muscle  needed  to  control  dispersion.  For  a  very  small  angle  of  attack  the 
control  moment  might  be  too  large  and  cause  large  roll  rates.  The  assumption, 
po  »  p+  must  be  upheld.  Also,  roll  rates  of  great  magnitudes  are  not 
practical . 

Three  types  of  control  are  examined:  derivative,  proportional, 
integral.  The  expressions  for  the  gains  are  shown  in  Table  1  for  the  case 
where  6o  31  iaQ.  The  magnitude  of  the  control  moment  is  the  product  of  the 
gain,  the  feedback  variable,  and  the  roll  moment  of  inertia. 

By  neglecting  the  higher  order  terms  in  Eq.  (30),  we  find  the  feedbacks 

to  be  independent  of  the  type  of  disturbance  (i.e.,  impulse,  step).  However, 

the  direction  of  the  disturbance  is  important  because  linear  superposition 

principles  do  not  hold.  For  an  arbitrary  moment  disturbance  with  components 

in  both  the  y  and  z  directions,  two  nonzero  feedbacks  are  required  to  satisfy 

Eq.  (34).  They  are  dependent  on  the  ratio  of  the  directional  components  of 

the  disturbance.  For  a  disturbance  about  one  axis,  only  a  single  feedback  is 

required.  The  feedback  variable  is  the  same  for  either  case,  m„  or  m,.  The 

y  * 

gain  depends  on  the  direction  of  the  disturbance  but  is  independent  of  its 
magnitude. 

An  initial  angle  of  attack  aQ  requires  the  nonzero  gain  to  be  B,  C,  or  E 
for  a  disturbance  either  in  the  y  or  z  direction.  If  the  disturbance  has 
components  in  both  directions,  A,  D,  or  F  must  also  be  nonzero.  These  results 
are  summarized  for  proportional  control  with  6  -  ia 


Case  I  -  mt  *  Oy 


P1  -  V 


P02[(o  -  l)2  -  4] 
2ao(o  +  1) 


Table  1.  Gains  for  5  -  ia 

o  o 


Control  Type 

Proportional 
P  *  bb+ 


Integral 

p  -  eJ<*+ 


Moment 

Direction 


Derivative 
•  • 
p  *  Ca 


Case  II  -  m*.  » 

t  z 

p2  -  b 2e 

-P02[(a  -  l)2  -  4] 
^2  aQ(a  +  l)(o-  1) 

Case  III  -  *  niy  +  imz 

P3  *  Pi^a  +  B33 


*3 


<B„  -  B.) 


I  +  Y 


2  v  2  1 


B.  + 


i  ,  21'  2  “2 

1  +  Y  1  +  y 


where 


Y  =*  -m  /m 

y  z 

If  the  initial  angle  of  attack  is  5q  »  Bq,  the  feedback  variable  in  Cases 
I  and  II  is  a.  The  single  feedback  has  always  zero  initial  value.  The 
perturbation  from  the  nonzero  Initial  angle  of  attack  does  not  have  to  be 
calculated. 

For  the  case  of  zero  initial  angle  of  attack,  -  0,  the  expression  for 
F(s),  Eq.  (27)  reduces  to 

F(s)  -  5+(s)  +15?  [S+(t)A(t)J 
-  \  2  [Mt)\2(t)]  + 


•  •  • 


( 


The  first-order  approximation  will  no  longer  hold  because  the  feedbacks  do  not 
appear  in  any  linear  terms.  A  second-order  approximation  of  F(s)  is  required 


F(s)  -  <5+(s)  +  i  5?  [6+(t)X(t)] 


(36) 


The  solution  for  zero  dispersion  is  found  from  Eq.  (31)  by  setting  F(-i)  *  0. 
The  first  term  in  Eq.  (36)  is  easily  evaluated.  With  s  *  -i  Eq.  (19)  becomes 


fi+(-i) 


a?(imt) 

Pq2(o  +  1) 


(37) 


The  higher  order  term  in  Eq.  (36)  must  be  evaluated  using  a  Laplace  trans¬ 
form  multiplication  theorem:5  "If  fj(t)  and  f2(t)  are  2? -transformable  func¬ 
tions  having  the  S?- transforms  Fj(s)  and  F2(s)  ,  respectively,  and  if 
Fj(s)  A  Aj(s)/B|(s)  is  a  rational  algebraic  fraction  having  q  first-order 
poles  and  no  others ,  then 


q  A  (s,  ) 

V  lf.(t)f2(t)J  -  l  F  (s  -  s  )  (38) 

k-i  vv 

By  substituting  f x ( t)  *  6( t)  and  f2(t)  -  X(t)  into  Eq.  (38)  and  setting  s  -  -i 
the  second-order  term  can  be  obtained. 

For  an  impulsive  moment  disturbance  the  solution  required  for  mc  -  my, 

®y  [2PqB  “  P02C  +  TTTI  +  1P04l4(0  +  1)_1)  ”  0  09) 


and  for  mt  *  im2 , 


n>z[2p0A  +  P02D  -  +  ip04l4<0  +  Dl  -  0  (40) 

No  solution  exists.  Addition  of  the  third-order  term  in  F(s)  gives  the  same 
results.  An  impulsive  moment  is  uncontrollable  by  this  approach.  There  is  no 
angle  of  attack  to  provide  the  lift  force  required  for  control. 
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For  a  step  moment  disturbance  about  the  y-  or  z-axis,  the  solution 
requires  one  nonzero  feedback.  The  expressions  for  the  gains  in  the  case  mt 
im2  are  shown  for  the  three  types  of  control. 

2.  PROPORTIONAL 


aPomz^/Qj  +  w/Q,  +  1/Q,) 


3 .  DERIVATIVE 


4 .  INTEGRAL 


where 


Po  (o  -  l)mz 


mz[o)/Q1  +  <d/Q2  +  (a  -  1 )/ Q3 ] 


4po4a5(2a)2  +  2poa)  -  po2) 


o+l 


4pQ4o5(2u)2  -  2pQu>  -  pQ2) 
o~+  1 


ai2po3[(o  -  l)2  -  4] 
a  +  1 


The  gains  are  inversely  proportional  to  the  magnitude  of  the 
disturbance.  The  controlling  force  is  the  angle-of-attack  lift  force  caused 
by  the  disturbance.  An  examination  of  Eqs.  (41),  (42),  and  (43)  shows  that 
the  gains  are  very  large.  With  the  values  o  »  24  and  pQ  *  lOx  rad/s  the  gains 
are  calculated  as  A  *  -4.78  x  10^,  D  -  1.64  x  104,  and  F  -  1.62  x  107.  Roll 
rates  with  these  feedbacks  are  abnormally  large  and  impractical.  The  initial 
assumption  of  |p0|  »  |p+|  no  longer  holds,  which  suggests  that  the  control  is 
ineffective. 


5 


NUMERICAL  EXAMPLES 


The  equations  of  motion,  Eqs.  (1)  and  (6),  can  be  expressed  by  equating 
the  real  and  imaginary  components 

a  +  va  +  p(2  -  y)8  +  [m2  -  p2(l  -  u)]a 

+[p+p(v-v)]3=m  (44 

m  y 

3  +  vp  -  p(2  -  p)a  +  [oj2  -  p2(l  -  u)]S 

-  [p  +  p(v  -  v  ))a  *  -m  (45 

m  z 

L0 

V  =  V  (0)  -  f  (8cos®  -  asin$)dt  (46 

y  y  tn  ' 

La 

V  «  V  (0)  -  J  (acos*  +  Bsin$)dt  (47 

z  z  m  * 

These  were  integrated  numerically  for  both  open-  and  closed-loop  responses  to 
impulse  and  step  disturbances  in  the  y-z  plane.  Cases  of  a  nonzero  initial 
angle  of  attack  are  shown.  The  moment  producing  the  initial  angle  of  attack 
is  found  from  the  equations  of  motion  in  steady  state. 

The  following  examples  show  undamped  behavior.  The  first-order  feedback 
gains  were  used  in  the  closed-loop  cases.  Inputs  and  system  parameters  are 
included  in  Table  2.  The  gains  are  summarized  in  Table  3. 

Open-  and  closed-loop  responses  to  an  impulsive  moment  about  the  y-axis 
are  shown  in  Figs.  1  and  2.  The  impulse  was  approximated  by  a  rectangular 
pulse  of  0.0001  s  duration.  The  magnitude  of  the  disturbance  was  such  that 
the  open-loop  oscillation  was  roughly  0.1-deg  amplitude.  The  closed-loop 
example  is  a  case  with  derivative  control  where  p  »  Co.  The  dispersion 
velocity  crossplots  are  shown  for  both  the  open-  and  closed-loop  in  Figs.  la 
and  2a,  respectively.  Fig.  lb  is  the  open-loop  response  of  o,  the  angle  of 
attack  component  directly  affected  by  the  nL,  disturbance.  The  closed-loop 


Table  2.  Inputs  and  System  Parameters 


m  =  22.68  kg 

L.  =  3.275  x  105  N/rad 
0 

pQ  =  10  rad/s 

a  “24 

u  =  0 

v  =0 

♦(0)  =  0 

v  (0)  =  -40.109  m/s 
w  (0)  *  0 
a(0)  -  5  deg 
8(0)  =  0  deg 

my  (steady  state)  =  1980.96  s  ~ 
mz  (steady  state)  =  0 


Table  3 .  Gains  for  «  5  deg 


A  B  C  D  E 


F 


my  .  0  118  ,755  330  0  -  324  ,410  0 

mz  0  10,325  -  3780  0  3,730,730  0 
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VELOCITY  Ivl,  m/s 


2a. 


Closed-Loop  Response  to  Moment  Impulse 
velocity  crossplot. 


behaviors  of  the  roll  rate  and  a  are  shown  in  Figs.  2b  and  2c.  The  open-loop 
mean  value  of  the  dispersion  velocity  w  is  -0.16  m/s.  With  the  closed-loop, 
the  roll  rate  oscillates  with  only  a  0.6  rad/s  amplitude,  and  the  mean  value 
of  the  dispersion  velocity  becomes  virtually  zero. 

Figures  3  and  4  show  the  response  to  a  step  moment  resulting  in  a  0.1-deg 
step  in  6.  The  open-loop  behavior  of  a  is  shown  in  Fig.  3b.  The  closed-loop 
example  demonstrates  proportional  control  of  the  form  p  3  Bg.  The  open  and 
closed-loop  crossplots  of  the  transverse  velocities  are  shown  in  Figs.  3a  and 
4a.  In  addition  the  behavior  of  the  roll  rate,  p,  and  the  feedback,  g,  are 
shown  in  Figs.  4b  and  4c  for  the  closed-loop  case.  The  mean  value  of  the 
dispersion  velocity  is  reduced  from  an  open-loop  value  of  0.8  m/s  to  a  closed- 
loop  value  of  zero.  The  roll  rate  amplitude  is  approximately  1.3  rad/s. 

C.  SINUSOIDAL  FEEDBACK 

The  missile  response  to  a  sinusoidal  roll  rate  modulation  of  the  form 

p  =  Pq  +  A  sin  qt  (48) 

is  examined  in  order  to  determine  the  parameters  A  and  q  required  for 
effective  dispersion  control.  It  is  convenient  for  this  analysis  to  use  the 
aeroballistic  equations  of  motion  for  the  complex  angle  of  attack,  Eq.  (2), 
which  with  the  roll  rate  behavior  of  Eq.  (48)  can  be  written 

5  +  (v  -  ipp)£  +  (u>2  -  ipv  )£ 

m 

A 

*  imt  exp  i[p0t  +  z(l  -  cos  qt)]  (49) 

A 

where  z  is  the  nondimens ional  ratio  A/q.  For  A  small  relative  to  pQ  the  roll 

coupling  on  the  left  side  of  Eq.  (49)  is  very  weak  because  p  and  v  are,  in 

m 

general,  small  terms.  For  a  first-order  approximation  to  the  roll  modulation 
we  assume  zero  damping  v  *  wm  *  0  and  p  -  0,  which  reduces  Eq.  (49)  to 

••  *S  A 

£  +  w  £  *  imt  exp  i[pQt  +  z(l  -  cos  qt)]  (50) 
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Fig.  4a.  Closed-Loop  Response  Co  Moment  Step.  Dispersion 
velocity  crossplot. 


With  these  approximations  it  is  apparent  that  the  response  to  an  impulsive 
moment  rat  with  zero  initial  angle  of  attack  is  not  controllable  because  the 
right  side  of  Eq.  (50)  is  defined  only  at  t  *  0  and  is  therefore  independent 
of  roll  rate.  Consider  the  response  to  a  moment  step  mt  *  -imz  “  ®  H(t).  Eq. 
(50)  can  be  written 

A 

2  -  i(z+P  t) 

5  +  u»  5  *  mH(t)e  [cos(zcos  qt) 

A 

-  i  sin(zcos  qt)]  (51) 


where  the  sinusoidal  terms  can  be  expanded  in  terms  of  Bessel  functions  with 
the  relations** 


cos(zcos  0)  *  J  (z)  +2  l  (-1)  J  (z)  cos(2k0) 
0  k*l  Zk 


jin( zcos  0)  -  2  l  (-1)  J_.  ,(z)  cos [(2k  +  1)0] 


Eq.  (51)  becomes 


2  —  i(z+p  t) 

?  +  (Vz)  ”  2J2(z)cos  2qt  +  2J,  (z)cos  4qt  -... 


A 

-i [2J j ( z)cos  qt  -  2J_( z)  cos  3  qt  +...]} 


_  i(z+pot) 

mH(t)e  U  (z)  +...+  2(-i)nJ  (z)cos  nqt  +  ...] 


The  Laplace  transform  of  Eq.  (54)  gives  for  £(s) 


?(s) 


—  iz  J  (z) 


me 


[— g-r  -  +... 

*•  c  —  in 


2  2  Ls  -  ip 

3  +  u>  ro 


2(-i)  J  ( z)( s  -  ip  ) 


(s  -  ipQ)2  +  (nq)2 


The  inverse  of  which  is 


c(t) 


aeizJ  (z)  ip  t 

O  r  0 


ip, 


[e  0  -  cos  wt - sin  wt  1  +...+  m  eizJ  (z)(-i)n 

u)  *  n 


2  2 
u>  -  p 


(  L  r  KP0+nq)t 

t  ■)  o  Le 


i(p  +nq) 


ui2  -  (pQ  +  nq)2 


-  cob  wt - - -  sin  cat]  + 

b)  J 


2  i2 

U)  -  (po  -  nq) 


i(po-nq)t 

l* 


-  COS  (lit  - 


i(p0  -  nq) 

- - -  sin  wt  J} 


The  transverse  velocity,  Eq.  (7),  is  obtained  by  integrating  £(t)  in  Eq. 
(56).  The  result  is 


L.«eizJ  U)  .  .  .  ipn  L  m  eizJ  (z)(-i) 

... .  „  8  o  re  sin  wt  o  .>  ,  .8  n 

V(t)  - - , - ~  [ - r- - - - 5-  (1  -  cos  wt)j  +...+  - - - 


,2  2.  1  iP 

■(w  -  p  )  ro 


u  2 

w 


1— 
w  -  (p  +  nq) 


i(po+nq)t  1(p  +  nq) 

re  -1  sin  wt  ro  /,  „\ 

~—2  t  -iCp-Vlq ? - w - 2 -  °  "  C08 

+  nq)  *o  ^  w 


l  re 


i(P0-nq)t 


i(p  -  nq) 

°  /  »  _  _ \ 


1  sin  ut 


l 


For  pQ  equal  to  integer  multiples  of  q,  a  singularity  occurs  at  n  3  ng  where 
p0  -  ngq  *  0.  If  we  let  pQ  -  ngq  3  e,  we  can  obtain  the  value  of  the 
singularity  as  e  +  0.  The  singular  term  in  Eq.  (57)  having  pQ  -  nq  in  the 
denominator  can  be  written 


lim 

e+o 


iet 
e _ 

ie 


lim 

e+o 


cos  et  +  isin  et 
ie 


1 


t 


(58 


We  are  interested  in  the  steady-state  or  d.c.  component  of  V(t),  which  is 
found  to  be 


^steady 

state 


tL89Te 

mp 


iz 


+  •  •  •+ 


2(-i)nJn(z) 

~  ,  2  2.  27 

1  -  (n  q  /pQ  ) 


n  a  . 

+  ...+  (-i)  SJn  (z)  (-i  -  itpQ)]  (59 

s 

_  2 

where  m  has  been  replaced  by  to  0t,  being  the  magnitude  of  the  trim  step. 

In  the  limit  as  A  and  z  0,  JQ(z)  1  and  Jn(z)  ♦  0,  and  Eq.  (59)  reduces  to 

the  open-loop  value  for  the  transverse  velocity  increment  due  to  a  trim  step. 

1 .  NUMERICAL  EXAMPLES 

Undamped  responses  to  a  step  moment  with  and  without  the  roll  modulation 

were  obtained  from  numerical  integration  of  the  equations  of  motion.  The 

open-loop  response  of  the  velocity  increments  is  shown  in  Fig.  5.  The  mean 
value,  calculated  from  Eq.  (59),  with  £  =  0  is  w  «  0.802  m/s.  The  two 
examples  with  roll  modulation  are  the  cases  where  q  3  pQ  (Fig.  6)  and  q  =  pQ/; 
(Fig.  7).  For  q  3  pQ  the  roll  control  is 

p  3  50x  cos  lOut  (60' 


The  amplitude  is  such  that  the  A  of  Eq.  (48)  is  roughly  equal  to  5.  The 
steady-state  response  of  the  transverse  velocity  components  is  predicted  from 
Eq.  (59)  with  n  3  0,1  and  n  3  (pn/q)  3  1. 
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Roll  Modulation  Response  to  Moment  Step 
w-component  of  dispersion  velocity. 
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Fig.  7b.  Roll  Modulation  Response  to  Moment  Step 
w-component  of  dispersion  velocity 


v  -  -0.317  t  +  0.094 


w  -  1.980  t  -  0.791 


For  q  *  p  /2  the  roll  control  is 


p  *  10 tt  cos  5itt 


Using  Eq.  (59)  and  taking  n  *  0,  1,  2  where  ng  ■  (pQ/q)  ■  2,  the  response  is 


calculated  to  be 


v  =■  0.458  t  -  0.090  (6 

w  -  0.193  t  -  0.865  (6 

The  mean  values  of  v  and  w  approximated  by  Eqs.  (61),  (62),  (64),  and 
(65)  agree  well  with  the  values  obtained  numerically  and  demonstrate  a 
capability  to  compensate  for  lift  nonaveraging  dispersion.  If  we  consider 
the  amplitude  A  of  the  roll  modulation  term  to  be  an  equivalent  control 
deflection,  we  can  write  a  control  law  of  the  form 

A  “  -Ar  -  Br  -  C  /  rdt  (6 

where  r  is  the  missile  cross-range  dispersion.  The  control  loop  is  shown  in 
Fig.  8.  A  first-order  approximation  to  the  feedback  gains  can  be  obtained 
from  the  steady-state  response  to  a  moment  step,  Eq.  (59).  For  small  values 
of  z  *  A/q,  the  Bessel  functions  Jn(z  )  for  n  1  are  small  relative  to 
JQ(z)  *  1  and  the  singularity  term  dominates  the  response.  For  the  case 
q  -  p0,  the  steady-state  transverse  velocity  is  approximately 


v  +  iw  »  - 


iLo0_e 
9  T 


[J  (z‘)  -  iJ.(t)  (j  -  ipct)] 


With  e 


*  1,  Jq(z)  *  1  and  J^(z)  “  z/2,  the  component  w  is 

w  »  -w  (1  -  tA/2)  (67) 

o 

where  wq  »  Lg0^/tnpo  is  the  open-loop  response.  If  we  let  r  »  z  and  consider 
proportional  control  only,  the  dispersion  z  from  Eqs.  (66)  and  (67)  is 
described  by 

w  Atz 

z  +  —  *  -v  (68) 

i  o 

The  solution  to  Eq.  (68)  is 

2  2 
Aw  t  „  Aw  t 

o  t  o 

z(t)  *  e  ^  [-w  /  e  ^  dt  +  z(o)]  (69) 

o 


which  with  the  change  of  variable  T  ”  t/Aw  /2  can  be  written 

o 


z(t) 


where 


(70) 


2  x  2 

D(t)  ■  e  /  eT  dt  (71) 

o 

The  function  D(t),  known  as  Dawson's  integral,  is  tabulated.**  We  can  select  a 
value  for  the  gain  A  to  limit  z(t)  to  some  prescribed  fraction  of  the  open- 
loop  dispersion  wQt.  Shown  in  Fig.  9  is  a  comparison  of  the  open-  and  closed- 
loop  response  for  a  gain  of  A  «  100,  which  from  Eq.  (70)  should  give  a  closed- 
loop  response  approximately  15%  of  the  open-loop  value  at  t  •  0.45  sec. 
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Fig.  9b.  Closed-Loop  Response  to  Moment  Step, 
z-component  of  dispersion. 


III.  CONCLUSIONS 


Transverse  velocity  resulting  from  moment  disturbances  that  perturb  the 
angle  of  attack  can  be  controlled  by  modulating  the  roll  rate  with  state 
variable  feedback.  Without  an  initial  trim  angle  of  attack,  the  dispersion 
velocity  is  uncontrollable.  Feedbacks  that  cause  the  transverse  velocity  to 
average  zero  can  be  calculated  by  linearizing  the  equations  of  motion  in  body- 
fixed  coordinates  and  obtaining  a  solution  in  the  Laplace  transform  domain 
without  the  need  to  invert  the  transformed  equations.  A  first-order  solution 
is  sufficient  and  shows  the  gains  to  be  independent  of  the  type  of  moment 
disturbance.  Numerical  examples  show  effective  control  with  only  small 
deviations  from  the  steady-state  roll  rate,  indicating  small  roll  moment 
requirements.  An  alternative  feedback  control  adjusts  the  amplitude  of  a 
small  harmonic  modulation  of  the  roll  rate  about  a  steady  value.  The  roll 
modulation  frequency  required  for  effective  control  is  derived  in  the 
aeroballistic  coordinate  system,  which  results  in  minimal  roll  coupling  with 
the  complex  angle-of-attack  motion.  A  quasisteady  solution  for  dispersion 
velocity  with  constant  amplitude  roll  modulation  is  then  used  with  a  slowly 
varying  amplitude  to  derive  a  control  law  to  limit  dispersion.  A  first-order 
approximation  for  the  feedback  gain  required  produces  effective  control  that 
agrees  well  with  predicted  values  for  the  magnitude  of  dispersion  relative  to 
the  open-loop  response  to  a  moment  step. 
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ABBREVIATIONS  AND  SYMBOLS 


A,B,C,D,E,F 

H(t) 

i 
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L 

L0 

in 

mt 

Am 
_  t 

m 

My,Mz 

P 

Po 

q 

r 


S 

t 

v 

V 

AV 


w 

z 

a 


feedback  gains 
unit  step  function 

AT 

pitch  or  yaw  moment  of  inertia 
roll  moment  of  inertia 
nth  order  Bessel  function 
roll  control  moment 
roll  moment  disturbance 
lift  force 

lift  force  derivative 
missile  mass 

complex  disturbance  moment,  my  +  im2 
magnitude  of  disturbance  moment 
moment  step 

pitch  and  yaw  disturbances 
roll  rate 
initial  roll  rate 
roll  modulation  rate 
y  +  iz 

Laplace  transform  variable 
aerodynamic  reference  area 
time 

y  component  of  transverse  velocity 
transverse  velocity  in  cross  plane,  v 
transverse  velocity  increment, 

Av  +  iAw 

z  component  of  transverse  velocity 
nondimensional  ratio  A/q;  coordinate 
angle  of  attack 
initial  angle  of  attack 


»  < 


.V.  A  • 

■a— JkdMh. 


+  iw 


angle-of-attack  perturbation 
angle  of  sideslip 
initial  angle  of  sideslip 
angle  of  sideslip  perturbation 

“Vmz 

complex  angle  of  attack,  6  +  ict 

initial  complex  angle  of  attack, 

8  +  ia 

o  o 

complex  angle  of  attack  perturbation 
unit  impulse  function 
complex  angle-of-attack  increment 
sinusoidal  control  parameter 
angle  of  attack  (Euler  angle) 

nondimens ional  roll  angle  pQt;  t/AwQ/2,  Eq.  (70) 

moment  of  inertia  ratio,  Ix/I 

aerodynamic  damping  parameter 

yaw  moment  damping  parameter 

complex  angle  of  attack,  8  +  ia 

»2/Po2 

roll  angle  relative  to  inertial  reference 
undamped  natural  pitch  frequency 


L ABO RATORY  OPERATIONS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  is  conducting 
experimental  and  theoretical  investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  space  systems.  Versatility 
and  flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  person¬ 
nel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  systems.  Expertise  in  the  latest  scientific  developments  is 
vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The  labora¬ 
tories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  environmental  hazards,  trace  detection;  spacecraft  structural 
mechanics,  contamination,  thermal  and  structural  control;  high  temperature 
thermomechanics,  gas  kinetics  and  radiation;  cw  and  pulsed  laser  development 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control, 
atmospheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmo¬ 
spheric  optics,  light  scattering,  state-specific  chemical  reactions  and  radia¬ 
tion  transport  in  rocket  plumes,  applied  laser  spectroscopy,  laser  chemistry, 
laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency 
standards,  and  environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence  and 
microelectronics  applications. 

Electronics  Research  Laboratory:  Microelectronics,  GaAs  low  noise  and 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quantum  electronics,  laser  communications,  Ildar,  and  electro- 
optics;  communication  sciences,  applied  electronics,  semiconductor  crystal  and 
device  physics,  radiometric  imaging;  millimeter  wave,  microwave  technology, 
and  RF  systems  research. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal 
matrix  composites,  polymers,  and  new  forms  of  carbon;  nondestructive  evalua¬ 
tion,  component  failure  analysis  and  reliability;  fracture  mechanics  and 
stress  corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and 
elevated  temperatures  as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  phys¬ 
ics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and 
ionospheric  physics,  density  and  composition  of  the  upper  atmosphere,  remote 
sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 


